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Abstract. This paper presents a generalization to image matching of the Hamiltonian ap- 
proach for planar curve matching developed in the context of group of diffeomorphisms. We 
propose an efficient framework to deal with discontinuous images in any dimension, for example 
2D or 3D. In this context, we give the structure of the initial momentum (which happens to 
be decomposed in a smooth part and a singular part) thanks to a derivation lemma interesting 
in itself. The second part develops a Hamiltonian interpretation of the variational problem, 
derived from the optimal control theory point of view. 
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1. Introduction 

This paper arose from the attempt to develop the multi-modal image matching in the framework 
of large deformation diffeomorphisms. Initiated by the work of Grenander, this context was deeply 
used since |Tro95) . especially with applications to computational anatomy. The method followed 
is the classical minization of an energy on the space of diffeomorphisms, which enables to compute 
geodesies on this space and to derive the evolution equations. In most of the papers, the group of 
diffeomorphisms acts on the support of the template; we add to this one a diffeomorphisms group 
action on the level set of the template. This action is a natural way to cope with the multi-modal 
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matching and could be, in a certain way, compared to the metamorphoses approach exposed in 
|TY05j : the metamorphoses are another way to act on the images but the goal is very different in 
our case. Matching in our context, is to find a couple (ry, (f)) which minimizes the energy 

Eir,,<j>)=D{Id,iv,^)f + ^\\r,oIoo^-^-It,,g\\l,, (1) 

with /q the initial function (or image) and Itarg the target function. Id is the identity map in the 
product of groups, cr is a calibration parameter. The distance D is obtained through a product of 
Riemannian metrics on the diffeomorphisms groups. 

All the complexity is then carried by the group of diffeomorphisms and its action: in the particular 
case of landmark matching, the geodesies are well described. The problem is reduced in this 
case to understand the geodesic fiow on a finite dimensional Riemannian manifold. It should be 
also emphasized that this problem can be seen as an optimal control problem. In [EMTYOH] . 
numerical implementation of gradient based methods are strongly developed through a semi- 
Lagrangian method for computing the geodesies. A Hamiltonian formulation can be adopted to 
provide efficient applications and computations through the use of the conservation of momenta. 
In [VMTYdi] . statistics are done on the initial momenta which is a relative signature of the 
target functions. The existence of geodesies from an initial momentum was deeply developed in 
|TY05j . but this work dealt only with smooth functions for /q (essentially H^) however with a 
very large class of momenta. An attempt to understand the structure of the momentum for an 
initial discontinuous function was done in the matching of planar curves in |GTL06) . 
We propose thereafter a framework to treat discontinuous functions in any dimension: the main 
point is to derive the energy function in this context. 

Finally, we chose to give a Hamiltonian interpretation of the equations which is the proper way 
to handle the conservation of momentum. This formulation includes the work done in jGTL06j 
but does not capture the landmark matching. The formulation we adopt gives a weak sense to 
the equations and we prove existence and uniqueness for the weak Hamiltonian equations within 
a large set of initial data. A word on the structure of initial data: the article on planar matching 
f |GTL06) ) focuses on Jordan curves. The main result is the existence for all time and uniqueness of 
Hamiltonian flow. The initial data are roughly a Jordan curve for the position variable and a vector 
fleld on this curve for the momentum. In our context, we choose four variables (/q, So,po, ^^o)- 
is the initial function with a set of dicontinuities Eq, Po is the momentum on the set Sq and Pq is 
the momentum for the smooth part of the initial function. This is a natural way to understand 
the problem and the choice to keep the set of discontinuity as a position variable can lead to larger 
applications than only choose /q as position variable. 

The paper is organized as follows. We start with a presentation of the framework underlying 
equation (U) . We present a key lemma concerning the data attachment term with respect to the ry 
and (j) variables. Its proof is postponed to the last section. Then we derive the geodesies equations 
and ensure the existence of a solution for all time from an initial momentum. In the second part of 
this paper, we give the weak formulation of the Hamiltonian equations, and deal with the existence 
and uniqueness for this Hamiltonian formulation. 
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2. Framework and notations 



2.1. The space of discontinuous images. Let n > 1 and M C M" a bounded open set 
difFeomorphic to the unit ball. 

We denote by BV{M) the set of functions of bounded variation. The reader is not supposed to 
have a broad knowledge of BV functions. Below, we restrict ourselves to a subset of BV functions 
which does not require the technical material of BV functions. However it is the most natural 
way to introduce our framework. Recall a definition of BV functions: 

Definition 1. A function / G L^{M) has bounded variation in M if 



In this case, Df is defined by Jj^j f div0 dx — — Jj^j Df (j) dx. 

Definition 2. We define Im(M) C BV{M) such that for each function f £ Im(M), there exists 
a partition of M in Lipschitz domains (J/i)ie[o.fc] for an integer k >0, and the restriction f\u^ is 



Remark 1. The extension theorem of Lipschitz function in M" enables us to consider that on 
each Ui, f^jj. is the restriction of a Lipschitz function defined on M". 

On the definition of a Lipschitz domain U : we use here (to shorten the previous definition) a large 
acceptation of Lipschitz domains which can be found in chapter 2 of [DZOlj . Namely, J7 is a hp 
domain if there exists a Lipschitz open set such that Q C U C 0,. In the proof of the derivation 
lemma [71 we give the classical definition of Lipschitz open set that we use above. In a nutshell, 
an open set is Lipschitz if for every point of the boundary there exists an affine basis of R" in 
which we can describe the boundary of the open set as the graph of a Lipschitz function on E"^^. 
We chose to deal with Lipschitz domains because it makes sense in the context of apphcation to 
images. 

Example 1. The most simple example is a piecewise constant function, f = 2_/i=i '^i-'-c/i ^^'^ 
a, e R. 

Remark 2. Our framework does not allow us to treat the discontinuities along a cusp, but we can 
deal with the corners respecting the Lipschitz condition. 

Let / £ Im(M), we denote by J/ the set of the jump part of /. As a BV{AI) function, we 
can write the distributional derivative of /: Df — V/ + D'^f + j(/)(a;)J{"^^L J/. V/ is the 
absolutely continuous part of the distributional derivative with respect to the Lebesgue measure 
and Dc is the Cantor part of the derivative. In other words, with the classical notations i{f){x) — 
if^ix) - f-{x))vf{x), where {f+,f-,Vf) : J/ R" x R" x 8"-^ is a Borel function. The 
functions /+ and are respectively defined as f^{x) — limtH^o+ f{x + t^f{x)) and f~{x) — 
limj^Q- + tvf{x)). Naturally, j{f) does not depend on the choice of the representation of Vf, 
in fact j{f) is homogeneous to the gradient. See for reference |Bra98j or [AFPOOj . In our case, 
the Cantor part is null from the definition. 
We then write for / G Im(M), 




Lipschitz . 
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Df = Wf+j{f){x):K'^-Kjf. (2) 

2.2. The space of deformations. We denote by Vm, <, >v a Hilbert space of square integrable 
vector fields on M, which can be continuously injected in {xq{M), ||.||p,oo), the vector space of 
with p > 1 vector fields which vanish on dM. Hence, there exists a constant cv such that for all 
vGV: 

||w||p,oo < Cy||w||y. 

Hence this Hilbert space is also a RKHS (Reproducing Kernel Hilbert Space), and we denote by 
kv{x, .)a the unique element of H which verifies for all v E H: < v[x),a >—< kv{x, .)a,v >h, 
where <, > is the euchdean scalar product and a a vector in M". This will enable an action on 
the support M. 

We denote by S", <, >s a Hilbert space of square integrable vector fields on R, as above. We denote 
by ks{x, .) its reproducing kernel. This will enable the action on the level set of the functions. 

Through the following paragraph, we recall the well-known properties on the fiow of such vector 
fields and its control. Most of them can be found in chapter 2 of [GlaOSj . and are elementary 
applications of Gronwall inequalities. (See Appendix B in |Eva98) ) 

Let V e L^([0, l],V), then with |Tro95j the flow is deflned: 

dt(blt = vtocbl,, (3) 
cl^Q = Id. (4) 

For all time t G [0, 1], 0o,t is a diffeomorphism of M and the application t rf^^o,* continuous 
and solution of the equation: 

dx(t)ot=Id+ / d^v ^^)Vs.dx4''"sds. (5) 

We dispose of the following controls, with respect to the vector fields; let u and v be two vector 
fields in G L^{[Q,l],V) and T < 1: 

llCt - '^atlloo < cy||w - u|lii[o,T] exp(cv|l'u|lii[o^T]), (6) 
\\d(t>Q,t - d(t)l^t\\oo < 2cv||w-u||ii[o,T]exp(cy||w||ii[o,T])- (7) 

And we have controls with respect to the time, with [s,t] C [0,r]: 

\\vr\\oodr<cv \\vr\\vdr, (8) 







< 


Uo.t- 




< 


Wd^o^t-d 


4>0,s\\oo 


< 



(9) 

\vr\\vdr. (10) 

with the constants C and C depending only on cy ■ Obviously these results are valid if S replaces 
V. In this case, we write 770,* for the fiow generated by St- With the group relation for the fiow, 

r]t,u ° ris,t = r]s,u- 
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The group we consider is the product group of all the diffeomorphisms we can obtain through the 
flow of u e L^{[0,1],V X S). 

We aim to minimize the following quantity, with /i the Lebesgue measure: 



3ivt,St) = - \\vt\\vdt+!- \\st\\i;dt+ |7?0.1%°'^oi(")-^targHrd/i(u), (H) 

^ Jo ^ JQ Jm 

with ryo,t the flow associated to sj. Remark that the metric we place on the product groups V x S 
is the product of the metric on each group which is represented by the flrst two terms in (flTI) . The 
functions Io,Itarg lie in Im(M). In one section below, we prove classically that there exists at least 
one solution and we derive the geodesic equations which give the form of the initial momentum. 

3. Derivation lemma 

This derivation lemma may be useful in many situations where discontinuities arise. Consider 
for exemple two Lipschitz open sets U and V. One may want to deform one of these open sets 
while the second remains unchanged (flgure below). The basic case is the following: 

Jt= [ xuo(j>;^dx^ fi{Vn(f>t{U)), 
Jv 

with ^ the Lebesgue measure. We answer to the differentiation of Jt, we obtain a sort of Stokes 
formula with a perturbation term. We discuss below a more general formula to apply in our 
context. The flnal result is the proposition [TJ 

Lemma 1. Let U,V two hounded Lipschitz domains of . Let X a Lipschitz vector field on 
M" and 4>t the associated flow. Finally, let g and f Lipschitz real functions on R". Consider the 
following quantity depending on t, 

Jt= f o (j)-^glvdii, 

Ut{u) 

where d^ is the Lebesgue measure, then 

dt\t=o+Jt = f - <Vf,X > glvd^i+ [ <X,n> fgiviX)d^ngu. (12) 

Ju JdU 

with lviX){y) = limeH^o+ '^viv + ^/ tf^-^ limit exists, elsewhere. And we denote by d^\gij 
the measure on dU and n the outer unit normal. 



As a corollary, we deduce: 
Corollary 1. We have, 

dt\t=o+KVr\MU))= I <X,n>iv{X)dfH9u, 

JdU 

with lviX){y) = limg^o+ '^viu + ^'^)! */ the limit exists, elsewhere. And we denote by d^^Qjj 
the measure on dU and n the outer unit normal. 

In this case, the derivation formula is a Stokes' formula in which one takes only into account the 
deformation viewed in V. 

Below is a flgure to illustrate the lemma: 
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Figure 1. Evolution of the area between two Lipschitz curves. (The arrows 
represent X along the boundary of U) 

Remark 3. We could generalize the lemma to finite intersection of Lipschitz domains, with the 
same scheme of the proof developed above. We gain hence generality which seems to he very natural 
for concrete applications. 

This generalization for Lipschitz domains is sufficient for the application we aim, and this appli- 
cation is presented in the paragraph below to derive the geodesic equations. Hence we present the 
corollary we use in the next paragraph. 

Theorem 1. Let {f,g) G Im{M)^ , X a Lipschitz vector field on K." and 4>t the associated flow. 

Jt= I f o (j)-^{x)g{x)dfi{x), 
Jm 

then the derivation of Jt is: 

dt\t=o+Jt = I -<V!,X>gdx~ /(/+ - < Uf.X > d'K^-\ (13) 

JM J 

with gx{x) := lim(^o+ g{4>t{x)) if the limit exists and if not, gx{x) = 0. 

Proof: Writing / as / = X^iLi f^xeUi where (C/i)j=i,- -,n is the partition in domains associated 
to /, and using the same expression for g, by linearity of integration, we fall in the case of the 
proposition [3l □ 

A last remark on the formulation of the lemma, we can rewrite the equation lfT3|) in a more compact 
form: 

dtit=o+Jt = - I {Df,X)g, 
Jm 

with Df the notation for the derivative for SBV function and g is the function defined above. 
Remark that fi a.e. g — g, these two functions differ on J/. 

4. Minimizing the energy functional 

The existence of geodesies is a classical fact, but in this framework the derivation of the geodesic 
equations did not appear to the author in the existing literature. With the metric introduced 
above, a geodesic in the product space is a product of geodesies. We chose to understand the 
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two geodesies separately for technical reasons. We could also have described the geodesies in the 
product space V x S, this point of view will be detailed in the Hamiltonian formulation of the 
equations. 

4.1. Existence and equations of geodesies. 

Theorem 2. Let (Iq, Itar) G I'm{M)^ , we consider the functional 3 on H ~ L'^{[0,l],VxS) defined 
in (fill) . There exists {v, s) ^ H such that s) = min(„_s)g// s). For such a minimizer, there 
exists {pa,Pb,Pc) e L2(Af,R") x L^{Ji„,W^) x L^{M,R) such that: 



(3st= Pciy)d[T]t.i]isi^y)ks{my),.)dn{y), (14) 

JM 

Xvt^ / kvi(l)o.tix), ■)[d(j)o^t]x^*iPa{x))dfi{x) 

JM 

+ / kv{(jJo,tix),.)[d(l)Qj]^^*{pb{x))dfj.\jj^{x), (15) 

with: 

It =m..i°io° 4>tfi, 

and Jig the jump set of Iq. More precisely for the (pa,Pc) we show, we have the equation: 

Pa{x) = -Jac((?!)o,i(a;))V|^/o Pc{x). (16) 

Proof: 

On the space H, the strong closed balls are compact for the weak topology. The functional 3 
is lower semi-continuous, so we obtain the existence of a minimizer. As reference for the weak 
topology |Bre94) . We find here |TY05j a proof that the flow is continuous for the weak topology, the 
main point to prove the semi-continuity: if (un,s„) ^ {u,v) in H, then ((/)o,i, ?7o,i) '-^ (0o,ii%,i)- 

We first differentiate w.r.t. the vector field s e L^([0,T], 5*), we denote by s a perturbation of s. 
Using the lemma flOl in appendix, we write: dsr]o,i{x) — jQ[dr]t,i]\rj„ ^(x)Sti'rio,tix))dt. We already 
introduce the kernel: 

[P < St, St >s + 2(/i - Itarg)[dr]t,i]\r,ot(i-(y)) < fcs(/|(y), .),«( >s dfj.{y)]dt = 0, 

Jo JM 

it leads to: 

Pst + [ 2(/i(y) - Itargiy))[dVtA]\my)ksimy), ^l^iv) = 0. 
JM 

With the notation pc — — 2(/i — Itarg), we have the first equation announced. 

For the second equation, we need the derivation lemma detailed in section [3l In order to use the 
lemma, we first need to develop the attachment term: 

byo.l % ° 0O.i(") - Itarg{u)\'^d^.{u) = / [Il f - 2I^Itarg + I?argd^J■, 
M JM 

Now, only the first two terms are involved in the derivation, and we apply the lemma to these two 
terms. (Actually the lemma is necessary only for the second term.) 
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Again, we have with the lemma flOl 

Vi := 9e(/)i,o(a;) = - / d{(j)t^o)4>i_t{x)iv[(l}i.tix)])dt. 



We consider the semi-derivation of (fTTj) at the minimum with respect to the displacement field v, 
we use the notations of SBV functions for the derivatives: 

A^' < vuvt >dt + J D[{r,)']Vi - 2 J Dl^VJl^g = 0, 
A / <vt,iJt>dt+ f 2 < V/i", Vi > {II - Itarg)dn + 

Jo JM 

OmY) - 2j(/r)) Vjtargd^''-^ = 0. 



As {IlY ^nd II have the same discontinuity set, the second integration is only over Ji^. 

We apply a version of the central lemma in |GTLQ6j which is detailed in appendix (see lemma [9]). 

g : L\[Q,l],V) ^ L\[QMV) L^M^W^) L^Ji^.W') 

We ensure that B := {(«, 2V ^{1^ ~ harg), {jm?) - 2j(/r)) harg), v G ^^([o^ i]^ y)} bounded. 
For each v, \Itarg\oo < \Itarg\oo- (This assumption could be weakened.) Whence we get with the 
lemma [9] the existence of Itarg G Conv(i?) (we observe that the Lebesgue part of Itarg is equal to 
Itarg, the modification is on the set J/j-), such that: 



A^ <Vt,Vt>dt+ J D[iI^ f]Vi - 2 J DllVjtarg = 0, 



A / < Vt.it >dt+ f 2V/l'(/r - Itarg)dtl + f j{[r,f)Vi - [ 2j (I^)ItargVi = 0. 

JO JM JJjV Jjrv 

Now, we aim to obtain the explicit geodesic equations by introducing the kernel, we denote by 
^i^) ■= / ^[(^r)^]^i ~2 J DliViItarg the pseudo derivative of the attachment term and we denote 
also: 

M^) iiihY) ° <l>l,o{x) - 2j(/o) O (j)lfi{x)itarg{x), 

which defines a normal vector field on Ji« . 

Mi) = / - / 2{h ~ Itarg)(y) < ky {4>l.t{v) , Wt^o)*^, ,iy){'^ ,iy)Io) .^t >V dfiiy) 
JO JM 

< A:y(?!)i,t(y),.)d(0t^o)^i.,(y)(A(y)), 1)4(2;) >v d^x\^^^ ^(j,^){y)dt. 

1 4>o,i(Jia) 

With the change of variable x = 0i,o(y), 

A{v) = I I 2(/i - /targ)(0o,i(a^))Jac(0oa(a;)) < fcv/((/)o,t(x), .)d((/)t,o)l /^)(V|^/o),{;t >y d/i 

Jo JM 

^tal^'f'^M < kv{<l>oA^),.)diM;o.i.)iMMx)),Vt >v d^^^j^Jt, (17) 
Jig \<J-(po,i[nx)\ 
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with a normal unit vector to J/^ in a; g J/^. Note that in the second term, the change of 

lac(0o i(a)) 
\d<i>o^i[n^)\ 



variable acts on the hypersurface Ji„ . This explains the term '^u^A,'^°fn^)\ which corresponds to the 



Jacobian term for the smooth part. 
We are done, 



Xvt = / 2A(0o,i(^))Jac(0o,i(a^))fcv(0o,t(x),.)d(fe);o,(.)(^l-^')^/^ 



M 



With pa{x) = 2A(0o,i(a^))Jac(0o,iW)V|.(/^) and pj,{x) = A(0o,i(a;))^i^7tSr' ^^^^ 
geodesic equations. □ 

These geodesic equations are a necessary condition for optimality. In the next paragraph, we show 
that if {pa,Pb,Pc) is given, we can reconstruct the geodesies. 

4.2. Reconstruction of geodesies with the initial momentum. We first demonstrate that 
if a vector field is a solution to the geodesic equations, then the norm is constant in time. 

Proposition 1. Constant speed curves in vector fields spaces 

If a vector field St is a solution of equation (fT4| and the kernel is differentiahle then ||st|p is 
constant. 

If a vector field vt is a solution of equation (flSl) and the kernel is differentiahle then ||wt|p is 
constant. 



Proof: 

We prove the first point: 

ll'St|ls= / / P{y')d[Vt,i]my')ks{my'),It{y))p{y)d[Vts]i^{y)dn{y')dfi{y). 

JM JM 

Remark that a.e. dt{d[r]t^i]i^=(y)) — —d[st]ii(y)d[rit^i]j!>(^yy This equation is obtained by a derivation 
of the group relation: d[77o,i]/|f — d[r]t^i o 770, t]/,^; and with the derivation of the equation (fT4)) : 

d[st]x^ / d[rits]^^yjPciy)diks{x,my))dy. 

J M 

As dst e i^([0, 1]) then the equation ^ proves that rf[?7t,i]/|(.) is absolutely continuous. As the 
space of absolutely continuous functions is an algebra, ||sf||| is also absolutely continuous. To 
obtain the result, it suffices to prove that the derivate vanishes a.e. 

dtWstf = - / / p{y')dstiiny'MvtA]^(y>Mitiy%^iy)))piy)d[vu^^^^ 

JM JM 

+ 11 p{y')d[i^tAiay')diks{It{y'))Jt{y))>tmy')^^^^ 

JM JM 

The second point is very similar. We underline that the equation (fT5|) is a particular case of 
the following, with a measure v which has a Lebesgue part and a singular part on the set J/q of 
discontinuities of the function /q. We also define: 

Pt{x) = {d[(t)o^t\l)~^{Pa{x)lx<fiJi^ +Pb{x)l^^.J,^^), 
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By the definition, 

||wt||^ y y Pt{x)kv{(t)Q,t{x),(t)o,t{y))pt{y)dv{x)dv{y). 
Remark that dtPt{x) = ^dvl o (j)Q t(x)pt{x), and we differentiate: 
9t||wt|P = - / j j pt{x)dkv{(j)oA{x),(t)aA{z))pt{z)kv{(j)o^t{x),(l)()^t{y)^^ 

Pt {x)dkv [4)0,1 [x) , (j}o,t {y))Pt {y)kv {4)0,1 (a;), 4)o,t {z))pt {z)dv{x)dv{y)dv{z) = 0. 

J J J 

□ 

This proposition is crucial to establish that the geodesies are defined for all time. Namely, we 
answer to existence and uniqueness of solutions to (the set J/^ but could be much more general 
than the discontinuity set of a function in Im(M)): 

?7o,t = Id+ I Suo riudu, 
Jo 

Pst{.) = / Pc{y)d[r]t,o\ii(y)ks{It{y)),-)d^i{y), 

JM 

4>o,t = Id+ 4>udu, 
Jo 

Xvt{.) = / kv{-,4>o,tix))[d4)oj]~^*{pa{x))d^i{x) 



M 

kvi;4o,tix))[dq)o,t]~^*iPbix))dii\j^^^ix)dt. (18) 

On purpose, this system of equations is decoupled in v and s. The proof of the next proposition 
treats both cases in the same time but it could be separated. 

Proposition 2. For T sufficiently small, the system of equations (fTSl) with 

{Pa,Pb,Pc) e Li(A/,R") X Li(J/„,M") X Li(M,M) 
has a unique solution if the kernel is differentiable and its first derivative is Lipschitz. 

Proof: 

We aim to apply the fixed point theorem on the Banach space i^([0, T],V x S). We estimate the 
Lipschitz coefficient of the following application: 

E: L'^{[Q,T],V X S) ^ L^i[0,T],V x S) 

{v,s) ^ (e("),^(s)), (19) 

with 



C(s)t = / Pc{x)d[r)t,o]i^(^)ksiIt{x)),.)d^i{x), 

J M 

^{v)t^ / kv{.,4o.t{x))[d(j)o,t]x^*{Pa{x))d^{x) 

JM 

kv{-, 4o,t{x))[d4)o,t]x^*{Pb{x))dfi\jj^ {x)dt. 

Jin 



HAMILTONIAN APPROACH TO GEODESIC IMAGE MATCHING 11 

For the space -^^^([0, T],V), if we have: 

mv)t - au)tf < M\\v - uhi[o,T], 

the result is then proven with Cauchy-Schwarz inequaUty: 



Uiv) - C(m)||l2[o,t] < vMT\\v ~ u\\l2[„^t]- 

This can be obtained with: 

\\^{v)t-au)tf - <av)t,av)t-au)t>-<au)t,av)t~au)t>, 

mv)t - iiu)t\\' < 2max(| < - > |, | < - > I)- 

For one of the two terms in the equation above: 

<av)t,av)t~a^)t> = J |[d0y-i*(p„(:r))[fc(0s,t(x),0;;,*(y))M0S,4"'*(pa(y)) 

- fc(0S,t(a:),(/.?f,i(y))M0?f,i]-^*(Pa(2/))]dA^(y)d/i(a;). 

On the unit ball of L'^{[0,T],V x S) denoted by B, and with the inequality we control the 
diffeomorphisms : 

ll<^o,t-<^S,tlU < cv\\v~u\\Li[o^T]eMcv), (20) 

\\d(j)^f~d(pgt\\oo < 2cy||w-u||Li[o,T]exp(cy). 

With the triangle inequality, we get: 

- Hcj^l^ix), (y))[d</>;^^,]-i*(p„(y))]|dA.(x)dA*(2/). 

On the unit ball B, we have: 

WltWo. < l + 2cyexp(cv), 
ll^o.i-^c'lloo <cv. 

Let Mk e K a bound for the kernel and its first derivative on the unit ball B. Such a constant 

exists thanks to the hypothesis on the kernel and its first derivative. 

A bound for the first term can be found with the second inequality of (|20| : 

2cv\\v - w||li[o,t] exp(cy)Mfc(l + 2cv exp(cy))||pa||||]36||, 

the second term is controlled with the first inequality of ((20| with the Lipschitz hypothesis on the 
kernel: 

cv\\v~ u\\l^[o^t] exp(cy)Mfe(l + 2cy exp(cy))||pa|| ||pf,|l. 

Finally we get, 

\\i{v)t - C(w)t|P < 6 cy ||u - u||li[o.t] exp(cy)Mfe(l + 2cy exp(cy))||pQ||||p6||. 

We have now concluded for the first component of the application S. For the second term, the 
proof is essentially the same, we do not give the details. □ 



12 



FRANCOIS-XAVIER VIALARD 



We have proven that there exists T > such that we have existence and uniqueness to the system 
(fTSl) . we prove now that the solutions are non-exploding i.e. we can choose T = +00 in the last 
proposition. This property shows that the associated riemannian manifold of infinite dimension is 
complete, since the exponential map is defined for all time. Without the hypothesis on the kernel, 
we can find simple counter-examples to this fact. 

Proposition 3. The solution proposition\M is defined for all time. 



Thanks to propostion [TJ we know that the norm of the solution ut is constant in time, which will 
enable the extension for all time. Consider a maximal solution with interval of definition [0, T] with 
T < +CG, then with the inequalities from ^ and after, we define the limit Ihnt^T (t>o,t '■= (f'o.T, 
since for all x, 0o,t(a^) is a Cauchy sequence. This is the same for limt^T d<po,tix). This limit is 
also a diffeomorphism, since we can define the limit of the inverse as well. The proof is the same 
to extend T]t for all time. 

We can then apply the proposition (existence for small time) to the current image It instead of 
Iq, we obtain diffeomorphisms 0o,s and 770, s in a neighborhood of 0, [0,e]. Composing with 0o,r 
and »7o,T, we extend the maximal solution on [0, T + e]. This is a contradiction. □ 

We decoupled the equations in s and v to give a simple proof of the existence in all time of the fiow. 
The formulation of ifTSj) implies the following formulation, which is the first step to understand 
the weak Hamiltonian formulation. If we have the system ifTSj) and the relation iflGj) . through the 
change of variable u — (/)o,t(a;), we get easily: 



5. Hamiltonian generalization 

In numerous papers on large deformation diffeomorphisms, the Hamiltonian framework arises. 
The simplest example is probably the Landmark matching problem for which the geodesic equa- 
tions and the Hamiltonian version of the evolution are well known ( |VMTY04) . | ATYOSj ) . Our 
goal is to provide an Hamiltonian interpretation of the initial variational problem. The main 
difference is that we want to write Hamiltonian equations in an infinite dimensional space, which 



Proof: 




(21) 



with. 



3a,c{(f>t,o)d[rit,o]mx)Po o 0t,o, 
-Pc{<l)o,i{x))d[rio^i]i^(^^-^3a,c{(l)Q^i{x)). 
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is roughly the space of images. The first step was done in |GTL06j where Hamiltonian equations 
were written on the representations of closed curves. Our work generalizes this approach to the 
space of images. We use the point of view of the optimal control theory (as it is developed in 
|ATY05j ) to formally introduce the Hamiltonian. Then, we state a weak Hamiltonian formulation 
of the equations obtained by the variational approach. We prove uniqueness for the solutions to 
these equations. At the end of this section, we discuss the existence of the solutions with the help 
of the existence of the solutions for the variational problem. 

In the whole section, we maintain our previous assumptions on the kernel for the existence of 
solutions for all time. 

5.1. Weak formulation. In this paragraph, we slightly modify the approach in order to develop 
the idea of decomposing an image in "more simple parts". Let introduce the position variables. 
We consider in the following that the discontinuity boundary is a position variable. Instead of 
considering the function It as the second position variable, we introduce a product space which 
can be projected on the space Im(M). Let (C/i, . . . , J7„) be a partition in Lipschitz domain of M. 
We denote by Sq — ^l^^i^Ui the union of the boundaries of the Lipschitz domains. We consider 
the projection: 

n 

p:Y[W^-'°°{M) Im(M) (22) 

n 

{P\=l.,...,n ^ / = (23) 

4=1 

Discontinuities give derivatives with a singular part. Maybe we could have treated this case 
adopting only the variable It, but we find the idea of decomposing an image into more simple 
parts rich enough to study the case. Observe that So is endowed with an important role in the 
definition of the projection: to write down a Hamiltonian system on the large space, we need to 
introduce the deformation of Eq and the deformation of each function in the product space. We 
will derive the Hamiltonian equations from this optimal control problem: (the position variable is 
Q and the control variable is U, c{U) is the instantaneous cost function) 

Q - {Q^)0<^<r = (S, {^)l<^<r) E ^^(Eo, M) X W^^^^MY, 

U = {v,s) eV xS, 

Q - f{Q,U)^{voQ'',{-{WQ\v)+s{Q'))i<,<r), 

c{U) = \\v\l + ^\s\l. 

The cotangent space D of the position variable contains F — L°°(E,K") x L^(M, R)'' C D. We 
write the formal minimized Hamiltonian of the control system on the subspace F, with P £ F: 

H{P,Q) = mmf (Po{x),Qo{x))dfi\j:,^{x) +^2 f P'Q^d^l{x) - c{U). (24) 

Minimizing in U , we obtain optimality conditions in (u, v) a minimizer such that for any pertur- 
bation {5v, Ss): 



14 



FRANCOIS-XAVIER VIALARD 



X{v,Sv) = 



JIM 



Using the kernel, it can be rewritten, 



Xv = f A;(Q°(rr;),.)P"dA*|E„(x)- V / kv{x, .)P'WQ'dpi{x) 



2A 



So So 

{y)VQ^ {y)ks (y, x)VQ' {x)F' (x)dii(x)dii{y) 



M JM 



- -2 J2 P\y)kv{Q\y),x)P\x)\/Q\x)dijL{x)diJi\^M 

+ h I I P'{y)ks{Q'{y),Q\x))p\x)di,{x)diJi{y). 

Now, we want to give a sense to the Hamiltonian equations. Vie [1, J"]: 



(25) 



Ps = y2 f ks{Q\.)P'dfi{x). (26) 
We deduce the expression of the Hamiltonian, 

HiP,Q) ^^[11 P"(x)fcy(QO(x),Q°(y))P"(y)dA.|s„(x)dA.|So(y) 



QO = dpoH{Pt,Qt){.), (27) 

Q\ = dp.H{PuQt)-iie[l,rl 

P° = -dQ.H{PuQt), 

PI = -dQ.H{Pt,Qt)^ie[\,r]. 

These derivatives should be understood as distributions, for * G C(|°° (Af, R) and u G C(f °° (M, R") 
and with the notation introduced in l|25p . Vie [1, f"]: 



9poi/(P,Q)(7.) = / (l'oQ>'(y),u(y))dAi|So(2/), (28) 

JT.0 

dp.HiP,Qm = [ ^iy){siQ\y))^{viy),VQ\y)))dfiiy), 

JM 

dQaH{P,Q){u) - / (H|QO(,)(«(y)),P"(y))dM|Eo(y), 



dQ,H{P,Qm = / vi/(y)[ds]|Q.(,)PXy)-(«(y),P'(2/)VvI/(y))d/i(y). 

JM 

Remark that only the last equation really needs to be defined as a distribution and not as a 
function. Now we can give a sense to the Hamiltonian equations but only in a weak sense: 
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Definition 3. An application x ^ C°([0, T], ^^(So, M) x VFi'°°(M)'- x ^^(Eo, M") x ^^(M, R)'-) is 
said to be a weak solution if it verifies for * e C^°°([0,T] x M,M.) and u e C^°°([0,T] x M,M"); 
(we denote x{t) = {Qt,Pt)-) 



(29) 
(30) 



/ f -dt^Qidfidt= [ dpiH{Pt,Qt){^)dt\/i€[l,r], 

Jo JM JO 

[ [ -a,.* Pidiidt = - ( dQ^HiPuQtW) dt^ie [l,r], 
Jo Jm Jo 

I I -dtuQ'ldii^^.dt^ [ dpoH{Pt,Qt)iu) dt, (31) 

Jo J'So Jo 

[ [ -dtuP°df,^^Jt=- I dQoH{Pt,Qt){u) dt. (32) 

Jo J'So JO 

5.2. Uniqueness of the weak solutions. In this paragraph, the uniqueness to the weak Hamil- 
tonian equations is proven, and the proof gives also the general form of the solutions. This form 
is closely related to the solution of the variational problem of the previous section. 

Theorem 3. Every weak solution is unique and there exists an element of L^{[0,T], V x S) which 
generates the flow (<?f)o,t j ??o,t) such that: 

Ql{x) = <^o,t(x),a:G So, (33) 
Qi{u) = r?o,toQo°</'t,o(w),w^^o,t(So),ie [0,n]., (34) 

and for the momentum variables: 

P^{x) = d[<t>o,t]-'*{P^{x)),xGi:o, (35) 
Pi{u) = Po*o0t,oJac((/)(,o)d[r?t,o]Q.(„),«^?!'o,t(So). (36) 

Proof: Let x ^ weak solution on [0,T], we introduce 

t - vt{-) = \ [ KQ'iix), .)P°dM|Eo(a^) - E / kv{x, .)PlVQ\d,x{x), 
Jt.0 Jm 

which Hes in L^([0,T], V). This vector field is uniquely determined by the weak solution. From 
the preHminaries, we deduce that (t)o,t{x) = j^Vs{(t)s,o{x))ds is well defined. We introduce also, 



t^st{.) = Y, I ks{Ql,.)Pidfi{x). 
■ 1 Jm 



For the same reasons, we can integrate the flow: rjo_t{x) = Sr{rirfi{x))dr is well deflned. In- 
troducing Ql{x) = rjt.o o Ql o (j)o^t for i G [1,^"], we obtain, with St o T]tfi{x) = dtr]tfi{x) and 
Vt o (po^tix) = dt(f>o,t{x): 
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'0 Jm Jo Jm 

/ -Vtfi ° Ql [dt'^] o (j)tfl^a.c{(j)tfi)diJbdt 

JM 
i-T 



/ / -%,o o Ql (5* [* o (t)tfl]- < V* o ^j_o, vt o (j)t^o >)i&c{(t)tfi)d^dt 
Jo Jm 



I (StiQi)- <VQi,Vt> +dr]tAQio(t>o,t))'^ diidt 
/o Jm 

= f f {St{Ql)-<VQl,Vt> + 
Jo Jm 

dvtA- < W o <f>o,t,vo 0o,t > +s{Q' o 0o,t)))* dfxdt. 

The cancelation of the equation above relies on the group relation of flows of vector flelds. We 
have the equality: 

St + dr]t,o{st o r]o,t) = 0, 
then the flrst and last terms cancel. The remaining terms cancel too because of the relations: 

VQj = #S,t(rfm,o(VQjo<^,,o)), 
vt + d(po,t{Vt o (j)tfi) = 0. 

Then, we conclude: 



Jo Jm 



Ql dndt = 0. 



Introducing ^{t,x) = X{t)j{x), with A G C+°°([0,T]) and 7 G C+°°{M), we have: 

lo -^'(*)(/m QidiJ,)dt = 0, hence: Jj^ 'y{x) Q\dn = Jj^ 7(2;) Q\dn, i.e. Q\ = Qi = Ql, and: 

Ql = m,toQo°4't,o- 

Now, we introduce for i G [l,r"], -P^ (•) = ^rf'rf ; this quantity is well defined because 

of the inversibility of the flow of St and Vt- Remark that '''^'^('^o.'(-)) — jg differentiable almost 

everywhere because Qq is Lipschitz on M . We want to prove that /o -St* PjM/xdt = with 
* G C^°°(M), which leads to: = -P^ = Pq'^ and then we are done. 

To prove the result, we flrst use the change of variable y = 4>o^t{x), this is a straightforward 
calculation, we will also use the equality: d[r]t,o]T^o ^(.)dr]o,t{-) = 1- 

r-T r pT 



I [ ^dtPl diidt = - [ [ Pidt^ o 0t,oJac(<^t,o) dtidt, 

Jo Jm Jo Jm 

[ [ ^dtPl dfidt = - / / — ^ dt^ o cj)t,o dfidt, 

Jo Jm Jo Jm HVtfilQii.) 



Jo Jm Jo Jm «W<,oJqj(.) «Wt,ojQj(.) 

+ P^^ o q^t^o dMvo,ty^o4,o.Mtxdt, (37) 
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The third term of the last equation can be rewritten: 

JM 

/ Pi-^o(j>t,o (< V(d[77o,t]Q-),-rf'/'t,o(ft) > +rf[st]Q.d[?7o,t]Q.o0o J d^^dt. (38) 

lO J M 

Now, we can apply the hypothesis on dP^, the first term of the expression is equal to: 

-Pidti-^^^) d^ldt = [ [ PI <v{^ocj,t.,od[m.t]Q^,o^oJ,vt> 

- * o 0t,o d[?]o,t]Qi^o4>oAst]QiPi dfidt, (39) 
Pi < V(^'o,/,t,od[r;o,t]Q.o0o J'^t > dfidt = f f Pi < d(/)*o(V«'|^,,J,z;t > d[T^o,t]Q'„o4>o t 

M ' JO JM 

+ PI < d05 ,V(d[77o,t]Q^),^^t > df^dt. (40) 
All the terms of the equation cancel together, so we obtain the result. 

With Pt{x) = d[(/)o,i]*P°(x) and = 0t,o ° Qo, we get the result for the last two terms of the 
system in the same way than the preceding equations, but it is even easier. □ 

Remark that we only have to suppose the weak solution is to obtain the result. More than 
uniqueness, we know that the weak solution "looks like" a variational solution of our initial prob- 
lem. 

5.3. On the existence of weak solutions. We consider in this section a Hamiltonian equation 
which includes our initial case, for which we have proved existence results. Namely, we can rewrite 
our result on the last section in terms of existence of weak solution of the system (|27|) . 

Proposition 4. Let {Ui, . . . , C/„) be a partition in Lipschitz domain of M , Sq = Uf^idUi. For 
any intial data, Iq G W^-°°{My, Qq = (Eo,/o) and Pq € L~(I]o,K") x L\M,Ry such that 
Supp{Pq) C Ui, then there exists a solution to the Hamiltonian equations. 

Remark that this solution has the same structure than a variational solution of our initial problem. 
In this case, all the momenta {i G [^,r]) can be viewed as one momentum: with Pt — J2i=i^tt 
we have all the information for the evolution of the system. 

Now, we can say a little bit more on the general Hamiltonian equations. We will not give a proof 
here of the existence if we relax the condition on the support of Pq, but the reader can convince 
himself that the existence is somehow a by-product of the last section. 

Summing up our work at this point, from a precise variational problem we obtain generalized 
Hamiltonian equations, for which we can prove results on existence and uniqueness. A natural 
question arises then, from what type of variational problems could appear these solutions? The 
answer could be based on the remark: because the decomposition we choose is the direct product 
of spaces, we can put a sort of product metric on it. A simple generalized minization is obtained 
by modifying the equation |(l]): 

E{r^, <f) = D{Id, {r^, <t>)f + ^ o (7^1^;.) o - iL-gWh. (41) 
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with for i € [l,r] 4 € W^'°°{M), and e Im(M). 

6. Conclusion 

The main point of this paper is the derivation lemma which may be of useful applications. 
This technical lemma gives a larger framework to develop the large deformation diffeomorphisms 
theory. The action on the level lines is far to be none of interest but we aim to obtain numerical 
implementations of the contrast term appHed to smooth images. Finally, the interpretation as a 
Hamiltonian system through optimal control theory ends up with giving a proper understanding 
of the momentum map. To go further, the technical lemma seems to be easily enlarged to rectifi- 
ables domains, and there may be a useful generalization to SBV functions. This would enable a 
generahzation of a part of this work to SBV functions. But to understand the weak Hamiltonian 
formulation would have been much more difficult within the SBV framework. From the numerical 
point of view, some algorithms that are currently developed to treat the evolution of curves could 
be used efficiently but they need strong developments. 

7. Proof of the lemma 
After recaUing some classical facts about Lipschitz functions, we prove the derivation lemma: 

Lemma 1. Let U,V two bounded Lipschitz domains o/R". Let X a Lipschitz vector field on 
K" and (j)t the associated flow. Finally, let g and f Lipschitz real functions on R". Consider the 
following quantity depending on t, 

Jt= fo (p-^glydfi, 

•l<t>t(U) 

where dfi is the Lebesgue measure, then 

dt\t=o+Jt = I -<Vf,X> glvdn+ I <X,n> fglv{X)dijL\au- (42) 

Ju JdU 

with lv{X){y) = lmi^^Q+ l-yiv + ^X), if the limit exists, elsewhere. And we denote by dfi^gu 
the measure on dU and n the outer unit normal of dU . 

We will use. 

Theorem 4. Rademacher's theorem 

Let f : U ^ M" a Lipschitz function defined on an open set U C K", then f is differentiable fi 
a.e. 

Theorem 5. Let F : 'R'^ x K" i-^ R'', Lipschitz continuous on a neighborhood of (vQ.yo) and 
F{vo,yo) = 0. Suppose that diF{vo,yo) exists and is invertible. Then there exists a neighborhood 
W of {vo,yo), on which there exists a function g : M" i-^ M.'^ such that in W : 

• aivo) = vq. 

• F{g{y),y) = for {g{y),y) G W. 

• laiy) - 9{yo)\ <c\y~yo\, 
withc=l + {Lip{F) + l)||9„F(^;o,2/o)"'||. 
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This theorem can be found in |PS03) and in a more general exposition than we will use hereafter. 
Now an obvious lemma of derivation under the integral, 

Lemma 2. Let f a Lipschitz function defined on an open set U C M". Let X a Lipschitz vector 
field on M" of compact support and 4>t the associated flow. Consider Jt = Ju f ° 4>tdx, then: 

dt\t=oJt= I <Vf,X>dx. (43) 
Ju 

Proof: Using Rademacher's theorem, this is a staightforward application of dominated convergence 
theorem. Note that under the condition that / is Lipschitz, if both / and V/ are integrable and 
X is a bounded vector field, we can relax the hypothesis of a compact support for the vector field, 
which is replaced here by the integrability condition. □ 

We will need the following characterization of derivation for real functions to prove the lemma |4] . 

Lemma 3. Let w : M" i-^ M a function, then w is difjerentiahle in x <E K" if and only if there 
exist f and g two functions and a neighborhood V of x, such that f{x) ~ g{x) and if y ^ V , 

9{y) < w{y) < fiy). (44) 

Proof: Suppose w differentiable in x, it suffices to prove that there exists / such that w < f 
in a neighborhood of x. (To obtain g, consider then —w.) we can suppose a: = 0, w{0) — and 
w'{0) = then limj^^^^: -^^j^ = 0. Hence there exists a continuous function v defined on -8(0, r > 0) 
such that: ^^j^ < v{y) and v{0) = 0. With the notation \x\ for the eucHdean norm in R", let 
TO(|a;|) — sup|j^|<|2.| v{y), we have m(|a:;|) > v{x) for x G B{0,r). The function m is non decreasing 
and continuous with m(0) = 0. At last, let f{x) = J^^^^ m{t)dt, then w{x) < |a;|m(|a;|) < f{x). 
Moreover, the fact / is is straightforward to verify. 

Suppose / and g are C^, and denote by /' and g' their derivative in x, then f'{x) = g'{x) since 
f ~ g ^ Q and has a minimum in x. On V we have: 

g{x + h) - g'{x).h < w{x + h) - g'{x).h < f{x + h) - f'{x).h. 

Hence, w{x + h) — g'{x).h = o{h) and w is differentiable. (Remark that we only use the fact that 
/ and g are differentiable in a;.) □ 

Using the lemma above, we study the deformation of an epigraph of a Lipschitz function under 
the action of a vector field, which leads to study the deformation of the graph of the function: 

Lemma 4. Let 4>t the flow of the vector field Lipschitz X on M" (with \\X\\ hounded on ]R"j. Let 
V = {{x,z) £ IR"~^ X Kj^ > w{x)} with w a Lipschitz function, and Wt{x) = inf{2;|(x, z) € (ptiV)}. 
Then, a.e. 

dt\t=awt{x) = - < Vw(x),pi{X{x,w{x))) > +p2{X{x,w{x))), 
with pi and p2 orthogonal projections respectively on M"^^ x and 0"^^ x M. 

Proof: Remark that Wt is well defined for all t by connexity reason, but Wt might be discontinuous 
for t large enough. However it is Lipschitz continuous for t in a neighborhood of 0: we first apply 
the implicit function theorem for Lipschitz maps to the function: 

F{x, t) = pi{(f>t{x, w{x))) - Xq, 
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note that dxF{xQ,0) = /c?|r„-i, so we obtain for each xq G M"^^ a function Xq : t t-^ xo{t) such 
that xq is Lipschitz and the equation F{x,t) = 4^ a; = XQ{t) on a neighborhood of (xo,0). Note 
that the impHcit function theorem in |PS03j gives only existence but not uniqueness. We develop 
now the uniqueness. The Lipschitz condition on w can be written with the cone property. Let 
{A, B) two points on the graph of w, then \yA — ys] < M\xa — xb], for a Lipschitz constant. This 
open condition is then verified in a neighborhood of {t = 0}. We see that, F{xi,t) = F{x2,t) 
implies 4>t{xi,w{xi)) = 4>t{x2,w{x2)), hence xi = x^. 

If w is C^, we get by implicit function theorem the first derivative of x^if): 

dt\t=oXt = ~pi{X{xt,w{xt))). 

We deduce, 

wtix) = p2i(f>t{xt,w{xt))), (45) 

and that 1 1— s- wt{xo) is a Lipschitz function for each xq. In the case, we get by differentiation 
of the equation (|45|) , 

dt^^^^wtixo) = - < Vw{x),pi{X{x,w{x))) > +p2{X{x,w{x))). 

We now observe that there is an obvious monotonicity in w of Wt. Indeed, if w < v then Wt < w*. 
We then use the lemma [3] to prove the result in the case w is Lipschitz . Let x such that w is 
upper and lower approximated by functions: let u and v such that u{x) = w{x) = v{x), and 
u < w < V. We obtain: 

^{ut{x) - u{x)) < j{wt{x) - w{x)) < j{vt{x) - v{x)), (46) 

We deduce the result: 

liin ^{wt{x) - w{x)) = - < \/w{x),pi{X{x, w{x))) > +p2{X{x, w{x))), 

for all the points of derivability of w, i.e. almost everywhere since w is Lipschitz. □ 

Now, we prove the following lemma, which can be seen as a consequence of the coarea formula. It 
will be used in the proposition [H 

Lemma 5. Let w : R" t-^ R, an a.e. differentiable function, A = w^^{{0}) and B = {x|Vw(x) 7^ 
0} then fi{AnB) = 0. 

Proof: For n = 1, the lemma is obvious because the point of AD B are isolated with Taylor 
formula. For n > 1, we generalize with Fubini's theorem: 

n 

AnB = [jAnBi, 

1=1 

with B, = {x e R"| < Vf{x),e, >=^ 0}. To prove the result, it suffices to see that n i?„) = 0. 
Consider (x,t) G R"^^ x R, fx{t) = f{x + 1) and Dx — x x R. Then we apply the case n = 1 to 
the function fx'- n{A n i?„ n Dx) = 0, and with Fubini's theorem, fi{A n Bn) = 0. □ 

Remark that this lemma can be appHed to a Lipschitz function. Below lies the fundamental step 
to prove the derivation lemma. 
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Proposition 5. Let w : R"^^ R a Lipschitz function. Let V :— {{x,y)\y > w{x)} and 
U := {{x,y)\y > 0}. Let X a Lipschitz vector field on M" and (pt the associated flow. Finally, 
let g and f Lipschitz real functions on R" of compact support. Consider the following quantity 
depending on t, 

Jt^ f o 4,-^glydfl, 

where d^i is the Lehesgue measure also denoted by dx, then 

dt\t=o+Jt^ -<'^f,X>glvdx+ < X,n > f giv{X)d^i\Qu ■ (47) 

JU JdU 

with lviX){y) = lime^o+ '^viv + ^/ tf^^ limit exists, elsewhere. And we denote by dn\Qu 
the measure on dU and n the outer unit normal. 




Figure 2. The main case 



Proof: The first case to treat is when w > 0, we can then integrate on V instead of 4>t{U). 

Jt^ f o (jy-^glyd^i, 
Jv 

we differentiate under the integral, we get: 

dt\t=QJt = / - < V/,X > glydfi, 

JU 

this is the formula because the second term is null. In the following, we have to use this case. 
To treat the general case, we first do a change of variable: 

p p+co 

Jt= I I fg o (j)t3ac{(l)t)dn. 

JR"-1 Jwt{x) + 

We introduce some notations: 

x^ = inax{Q, x) ^ H{x), 
VH{x){v) = 0,x<0, or x^Oetv <0, 
VH{x){v) — V , elsewhere, 

wtix) = M{z\{x,z) &<t,^t{V)}. 



Let pi and p2 the orthogonal projections respectively on M" ^ x {0} and 0" 
lemma [31 

dt\t=o+Wt{x) =< Vw{x),pi{X{x,w{x))) > ~p2{X{x,w{x))). 



With the 
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Consequently, 

dt\t=a+wt{x)+ ^ VH{w{x)){< Vw(a;),pi(X(x, w(a;))) > -p2iX{x,wix)))). 
Using fiy)i< Vg,X > +div(X)g) = div(/gX) ~g<Vf,X>, we get: 

dt\t=a+Jt^ / ~dt\t=()+wt{x)+f{x,w{x)+)g{x,w{x)+)dx 

< V f{x, z), X(x, z) > g{x, z)dxdz 

wt(x)+ 

+ fg<X,n>dfj.\9(^ijrtV)- 
Jd{unv) 

Here n is the outer unit normal of d{U n V). Rewrite the last term: 

/ fg <X,7i> dfj,\d!^unv) = / fg <X,n> dfj,\ov 

Jdiunv) Jdvnu 

f{x, 0)g{x, 0) < X,n> dx. 

u)-i(]-oo,0]) 

In a neighborhood C oi x such that w{x) > 0, we have demonstrated that the formula holds, so 
the first term in the equation above is equal to: 



fg<X,ji>dfj.\9v- / dt\t=omixy f{x,0)g{x,0)dx. 

IdVnU Jto-i(]-cx3,0]) 

Moreover, on the set F — {x : w{x) = 0}, we have, with the lemmaEJ a.e. Vw = 0. Then, we have: 
dt\t=o+'^t{x)'^ =< X,n >+ , a.e. On the set G = {x : w{x) < 0}, we have: dt^t^()+wt{x)~^ = 0. 
We now get the result with: 

{< X,n> -dt\t=owt{x)~^)f{x,0)g{x,0)dx = 

u;~i(]-oo,0]) 



/ 



lv{X){x,0) <X,n> f{x,0)g{x,0)dx. 

/u)-i(]-oo,0]) 

Indeed, ii < X,n >y^ the result is straightforward because the limit exists in the definition of 
ly. If < X,n >— 0, the contribution is null. □ 

Our goal is to prove the formula for Lipschitz open sets, we present some definitions. 
Definition 4. 

An open set U ^9 o/M" is said to be locally Lipschitz if for each x e dU , there exist: 

• an affine isometry I, ofW", 

• an open neighborhood V{x) of x, 

• a Lipschitz function w defined on M"^^ with a Lipschitz constant K[x) 
such that, 

I{V{x) n f/) = I{V{x)) n {{x,y) e M""^ x M\y > w{x)} 

If the constant K{x) can be chosen independent of x, U is said to be Lipschitz. 
Remark 4. 

(1) An open bounded set o/K." which is locally Lipschitz is also Lipschitz . 

(2) By Rademacher's theorem, the outer unit normal n{x) exists for a.e. x G dU . 
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(3) We will say that (V{x),I) trivializes the Lipschitz domain in x. 

The three lemmas below prove that one can describe a Lipschitz domain in many systems of 
coordinates. This is a key point to understand the two boundaries at a point of intersection and 
enables to use the proposition O 

Lemma 6. Let ij^ o, diffeomorphism o/M" and V a Lipschitz domain, then V-'Cl^) is a Lipschitz 
domain. 

Proof: The proof is straightforward with the characterization of Lipschitz domains with the uni- 
form cone property, which can be found in |DZ01) . □ 

Lemma 7. Let (ei, . . . ,e„) an orthonormal basis o/M", and w a Lipschitz function defined on 
]^n-i of Lipschitz constant M. Let U the Lipschitz open set which is above the graph of w: 

U ■■= {{x,y) e M""i X M|y > w{x)}, 

then for each n in the open cone C := {n = {x,y) £ M"^^ x R+|?/ > A/|a;|} one can trivialize 
the boundary of U through the graph of a function defined on (with an orthonormal basis). 
Moreover, this function is Lipschitz . 




Figure 3. Trivializing with respect to the orthogonal hyperplan to n. 

Proof: Note that, n can be represented with the angle between the hyperplan and the vector 
orthogonal to this hyperplan, and also the Lipschitz constant can be represented as the tangente 
of such an angle. Let two points a,b £ R" which belong to the dU , then b — a and n are not 
coHnear, because of the Lipschitz property. As a consequence, dU is defined as the graph of a 
function w on . And, one can verify that, if n is normalized, a Lipschitz constant for w is equal 
to: tan(\9i — 92\), if Oi is the angle of n and 62 associated to the Lipschitz constant. □ 

Lemma 8. Let U C M" a Lipschitz domain with G dU , then there exist a neighborhood V ofO, 
w a Lipschitz function defined on R"~^ and a linear transformation A such that: (ei, . . . , e„_i) C 
Ker{A-Id) and 

V n A{U) ^Vr\ {{x, y) e R" X M\y > w{x)}. 
Here is the illustration of the idea driving the proof. 

Proof: In some coordinates, we write U as the epigraph of a Lipschitz function w on a hyperplan 
if in a neighborhood W of 0. We note by n a normal vector to H. We face two cases: 
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A(M) 



Figure 4. The linear transformation A 



• If < n, e„ >7^ 0, let A G L(R") defined by: A{n) = e„ and (ei, . . . , e„_i) = Ker{A-Id). Denote 
by p2 the orthogonal projection on (ei, . . . , e„_i). Let M £ G defined by its projections z on H: 
M — z + w{z)n, by definition of A, A{M) ~ {p2{z),w{z)) + A{z — P2{z))- But z — P2{z) = X{z)n 
with A € (M")' so we obtain: A{M) = {p2{z),w{z) + A). Also P2\H is a linear isomorphism, we 
note the inverse p^^, then with the change of variable x = P2{z), we get 



w o P2 is clearly Lipschitz and we obtain the lemma in this case. 

• If < n, e,i 0, we can choose by lemma [7] another system of coordinates for which we fall in 
the first case, and the lemma is demonstrated. □ 

We present a smooth (C^) version of the derivation lemma. 

Proposition 6. Let U a bounded domain o/R" and V a bounded Lipschitz domain. Let X a 
Lipschitz vector field on R" and 4>t the associated flow. Finally, let g and f Lipschitz real functions 
on R" . Consider the following quantity depending on t, 



with iviX){y) = \im^^Qly{y + eX), if the limit exists, elsewhere. And we denote by dfi^Qu the 
measure on dU and n the outer unit normal. 

Proof: Let K a Lipschitz constant of V . Applying the definition of a domain with the compact 
boundary of J7, there exist a finite covering Wi, . . . , Wk of dU with open balls and . . . , VPfc) 
diffeomorphisms such that for each i e [1, A:], 



G = {(x, X{p2\x)) + wop-\x))\x e P2{H n W)}. 




(48) 



^i{W, nU) = {{x, y) e M"-' X M\y > 0} n ^^{W,) 



also one has: 



^^{w^ n du) = M""^ X {0} n ^'i(Wi). 

Let {9o, . . . ,9n, 9k+i) a partition of unity associated to the family 

(Wo = M" \ C7, VFi, . . . , Wk,Wk+i = U). 



It means: 
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• < 0, < 1, Vi e [0, fc + 1] and E»6[o,fe+i] O^^lon W\ 

• Supp 9i C Wi for ie [1, fc + 1]. 

• Supp Oq C R" \ U. 

Through the change of variable y = 4>t{x), the quantity is: 

fc+i » 

Jt ■^^'^ ° ° <j)tJa.c{(j)t)dfi. 

Four cases appear: 

• c7 n c 

• VK, n[/cM"\y 

• Y r^W^(lU 

In the first case, the formula is the result of the derivation under the integral, which is allowed 
because g is Lipschitz. 

9t|t=o+^t= / /(j/)(< V.g,X>+div(X).g)dAt, 

with, f{y){< V.g, X > +div(X)g) = div(/gX) — g < V/, X >, and applying Stokes theorem true 

for a rectifiable open set and Lipschitz functions, we obtain the result. 

In the second case, the quantity is null for t sufficiently small. So the formula is obvious. 

In the third case, we can integrate on V instead of 4>tiU n Wi): 

Jt^ f o (jj-^glydfi, 
Jv 

we differentiate under the integral, we get: 

dt\t=o+Jt^ / - <Vf,X > glydfi, 
JU 

because the second term of the formula is null. 

We deal hereafter with the last case: as is a diffeomorphism, '^i{V) is also Lipschitz. 
Consequently, we can find a finite covering Bi, . . . , B„i of Wi , for which one of the following 
conditions holds: 

• B,CV 

• B,cM."\V 

• BiH dV ^ and there exists / such as /) trivializes the Lipschitz domain ^i(V^). 

In the first two cases, we have already demonstrated that the formula is true. 
With the lemmalHl we know that after a linear transformation A which is the identity on M"^^ x {0}, 
the Lipschitz domain can be represented as the epigraph of a Lipschitz function defined on E"^^. 
We replace "^ihyAc^i — We then have the following situations: 

^iw,nv)^{iz,t)\t>w{z)}, 



or 



^{w,nv)^{{z,t)\t<w{z)}, 
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with w : R"^^ n -6(0, p) i— > R a Lipschitz function. This situation (or the symetric situation which 
is essentially the same) is treated in the proposition [H □ 

We generalize the proposition [6] to the case of Lipschitz domains, we need some additional results 
of approximation: 

Theorem 6. approximation 

Let f : M" i— > R a Lipschitz function. Then for each e > 0, there exists a function f : R" i-^ R 
such that: 

/i({x|/» ^ f{x)orDf{x) ^ Df{x)}) < e. 

In addition, 

sup\Df\<CLip{f), 

R" 

for some constant C depending only on n. 
See the proof of |EG92j . 

Remark 5. A direct consequence of the theorem is that we have, 

11/ -/loo < 2m£ix{C,l)Lip{f)V^^e^ . 

On each cube of volume e there exists a point where the two functions are equal, then we deduce 
easily the claimed bound. Thus, we get also 

l^i{ix,y)\f{x) <y< fix) or f{x) <y< f{x)}) < 2max(C, 1) Lip{f)^/n - le^+^ . 
We deduce a corollary: 

Corollary 2. Let U a bounded Lipschitz domain, for each e > there exists V a domain such 
that, S = U \ V U V \ U , is a rectifiable open set verifying: 

KS) < e (49) 
Jf"^i(a5)<e. (50) 

Proof: We just present the main points for the proof of the corollary. 

By compacity of dU, there exists a finite open covering {Vi, . . . , Vk) of dU, such that for each open 
set we can trivialize the boundary. On each open set Vi , we have by previous theorem a Lipschitz 
application gi : dU nVi i—* R" which gives a hypersurface. We have 

D<"-\{x £dUnV,\x^ gix)}) < e. 

Moreover we can assume that this covering satisfies the following property. Let < 77 < e and 
Z := U^^jVi n Vj, JC"~i(9[/ n Z) < 77. We thus obtain an application g : dU ^ R" which is 
Lipschitz (dU is endowed with the induced metric by the euclidean metric on R") and is the 
boundary of a domain V, for which we have: 

3<.'^-\{x e dU\x ^ g{x)}) < (fc + l)e. 

Then, J{"^^{dS) < 2 Lip{g)"^^{k + l)e. And also with the same argument given in the pre- 
ceding remark, there exists a constant K such that, fi{S) < K Lip(g) {{k + l)e) "-1 , with 
K = y/n^Up{dU). □ 
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We now turn to the proof of the lemma [TJ 

Proof: We use the corollary [21 let a domain for e as in the corollary. Let Mi a constant 
such that in a compact neighborhood of U, \ f ° (f't — f\ < Mit, g < M2, |/| < M3 and \X\ < K. 
We have, with 

S,^ ^{U,,U) = U,\U\JU\U„ 
We denote by 6* = — ly, , so we have (triangular inequality for the second inequation): 



\{Jt{U) - MU)) ~ {Jt{U,) ~ UU,))\ < / \eo(j)~^foq^-^g~0fg\dtl, 



V 

\{Jt{U) - UU)) - {Jt{U,) - MU,))\ < f |0o0-l(/o0-lg-/g)|d^ + 

Jv 

\{eo4>^' ^e)fg\di,, 

V 

liMu) - MU)) - (Jtiu,) ~ Mu,))\ < [ |/-/o0-i||5|iydM 

Jvr\4>t{s,) 

JvnA{MSe),S,) 

\{Jt{u) - Mu)) - {Jtiu,) - Mu,))\ < tAhM2fi{MSe) n v) + M3M2fi{vnA{MS.),s,)) 

We first treat the last term. We claim that, for sq > such that Lip((/)t) < 2, we have, for 
t e [-So, So], 

ti{A{MSe),S,)) < imax(2,M)"?{"-i(a5,))- 
Introduce ^ : {t, x) G [-so, sq] x R" ^ (j>t{x) e R". 

We have Lip(«') < max(2,M), and ?f"([0,i] x dS^)) = m''-'^{dSe). Hence, 3i:"(^'([-so, sq] x 
dSe)) < tmax(2, Af)"JC""i((9S'e)). To finish, we prove that: 

s<t 

Let ze A(0t(5,),5e), 

• Suppose z ^ Se, there exists x £ Sc such that (t>t{x) — z. The map c : s e [0,t] 1-^ (j)s{x) 
verifies c(0) = x G S and c(t) — z ^ S. By connexity, there exists u e [0,t], such that 
c{u) e dS^. By composition of fiow, 4>t^uic{u)) = z. 

• Suppose z ^ 0t(S'e), there exists x £ ^((S'e) such that <l)-t{x) ~ z. The map c : s e [0,<] 1-^ 
(j)-^s{x) verifies c(0) € 4't{Se) and c(i) = z ^ (f>t{Se). By connexity, there exists w € [0, i], 
such that c(u) G d(j)t{Se). By composition of fiow, to = e dS^ and obviously. 

We give a bound for the first term in the same neighborhood for t e [— sq, so], 

^l{MS,) nV)< fiiMSe)) < 2"tefi{S,). 

Consequently, 

lim sup \l[{Jt{U) - MU) - (JtiU,) - MUM < MiM2fi{S,) + Ms^^'-^dS,)). 
We can now obtain the conclusion. Let e > 0, 
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lim sup \-[{Jt{U) - MU) - {Jt{U,) - MUM < (A/iAf2 + Ah)e. 



We use now the formula already demonstrated for domains, 



lim sup lliiMU) ~Jo{U)- [ ~<Vf,X> glvdfi + 




f <X,n> fgiv{X)dn\ouM < (MiAh + Ms)^, 



JdUe 



and the result is proven. □ 



8. Appendix 



8.1. Central lemma of |GTL06) . We present here a different version of the lemma, which is 
essentially the same, but from another point of view. 

Lemma 9. Let H a Hilbert space and B a non-empty bounded subset of E a Hilbert space such that 
there exists a continuous linear application g : H ^ E. Assume that for any a <E H , there exists 
ha E B such that (ba,g{a)) > 0. Then, there exists b e Conv(_B) such that (b,g{a)) = 0, Va G H . 

Proof: We denote by Ho = g{H). Let p the orthogonal projection on ifo and Z = Conv(i?) is a 
non-empty closed bounded convex subset of H, whence weakly compact. As p is weakly continuous 
and linear, p{Z) is a weakly compact convex subset and thus strongly closed. From the projection 
theorem on closed convex subset, there exists b e piZ) such that: |6| = iidcez \c\ and {b,c—b) > 
for c £ p{Z). As a direct consequence, we have also: 



The element b lies in the adherence of g{H) then there exists a sequence &„ G g{H) such that 
lim 5,1 = b. From the hypothesis, there exists u„ € B such that: {un,—bn) > 0. As i? is bounded, 
lim(u„,6— bn) = 0, hence limsup(u„, < 0. By ((5T1) . we get limsup(w„,6) > (&, 6). As a result, 
(6, 6) < and 6 = 0. By definition, there exists w e Z, p{v) — b — 0. Now, v e Hq and 4>{.) — {v, .) 
gives the result. □ 

8.2. A short lemma. We give here a short proof of the perturbation of the flow of a time 
dependent vector field with respect to the vector field. We assume in the proposition that the 
involved vector fields are but it can be proven with weaker assumptions on the regularity of 
vector fields. (See [GlaOSj . for a detailed proof following another method.) 

Lemma 10. Let ut and vt be two time dependent vector fields on R", and denote by (I)q j the 
flow generated by the vector field ut + evt, then we have: 

de<j)o,iix) = / [d<j)t,i]4>o,tix)'"i(t^o,t{x))dt. 



Jo 

Proof: Introduce the notation At G M" defined by: At{(po,t{x)) = df_(t>Q t{x). Deriving this expres- 
sion with respect to the time variable: 

4:A(</'o,t(a;)) = dut{At{(l)Q^t{x))) +Vt{(j)Q^t{x)), 



{b,u-b)>0 Vit e Z. 



(51) 
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Remark that the expression above can be written as (with L the Lie derivative): 

By integration in time, we obtain the result. □ 
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